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Background: Neutrino-nucleus quasi-elastic scattering is crucial to interpret the neutrino oscilla- 
tion results in long baseline neutrino experiments. There are rather large uncertainties in the cross 
section, due to insufficient knowledge on the role of two-body weak currents. Purpose: Determine 
the role of two-body weak currents in neutrino-deuteron quasi-elastic scattering up to GeV ener- 
gies. Methods: Calculate cross sections for inclusive neutrino scattering off deuteron induced by 
neutral and charge-changing weak currents, from threshold up to GeV energies, using the Argonne 
«i8 potential and consistent nuclear electroweak currents with one- and two-body terms. Results: 
Two-body contributions are found to be small, and increase the cross sections obtained with one- 
body currents by less than 10% over the whole range of energies. Total cross sections obtained by 
describing the final two-nucleon states with plane waves differ negligibly, for neutrino energies > 500 
MeV, from those in which interaction effects in these states are fully accounted for. The sensitivity 
of the calculated cross sections to different models for the two-nucleon potential and/or two-body 
terms in the weak current is found to be weak. Comparing cross sections to those obtained in a 
naive model in which the deuteron is taken to consist of a free proton and neutron at rest, nuclear 
structure effects are illustrated to be non-negligible. Conclusion: Contributions of two-body cur- 
rents in neutrino-deuteron quasi-elastic scattering up to GeV are found to be smaller than 10%. 
Finally, it should be stressed that the results reported in this work do not include pion production 
channels. 

PACS numbers: 25.10.+S, 25.30.Pt 

I. INTRODUCTION 

In last few years, inclusive neutrino scattering from nuclear targets has become a hot topic. Interest has been spurred 
by the anomaly observed in recent neutrino quasi-elastic scattering data on 12 C [1, 2], i.e. the excess, at relatively 
low energy, of measured cross section relative to theoretical calculations. Analyses based on these calculations have 
led to speculations that our present understanding of the nuclear response to charge-changing weak probes may be 
incomplete [3], and, in particular, that the momentum transfer dependence of the axial form factor of the nucleon, 
specifically the cutoff value of its dipole parameterization [4] , may be quite different from that obtained from analyses 
of pion electro-production data [5] and measurements of the reaction n{v fl ,^)p in the deuteron at quasi-elastic 
kinematics [6, 7] and of v^p and V^p elastic scattering [8] (A^ ~ 1.20 GeV versus A^ ~ 1 GeV). However, it should 
be emphasized that the calculations on which these analyses are based use rather crude models of nuclear structure — 
Fermi gas or local density approximations of the nuclear matter spectral function — as well as simplistic treatments of 
the reaction mechanism, and should therefore be viewed with skepticism. 

In this paper, we calculate cross sections for inclusive neutrino scattering off deuteron in a wide energy range, from 
threshold up to 1 GeV. The motivations for undertaking such a work are twofold. The first is to provide a benchmark 
for studies of electro-weak inclusive response in light nuclei we intend to carry out in the near future. The second 
motivation has to do with plans [9], still under development, to determine the neutrino flux in accelerator-based 
experiments from measurements of inclusive cross sections on the deuteron. In particular, in charged-current neutrino 
capture on deuteron, the final states ppl~ can be measured, in principle, very well. Clearly, accurate predictions for 
these cross sections are crucial for a reliable determination of the flux. 

A number of studies of neutrino-deuteron scattering at low and intermediate energies (< 150 MeV) were carried 
out in the past decades, see Rcf. [10] for a review of work done up to the mid 1990's. These efforts culminated in the 
Nakamura et al.'s 2001 and 2002 calculations of the cross sections for neutrino disintegration of the deuteron induced 
by neutral and charge-changing weak currents. These calculations were based on bound- and scattering-state wave 
functions obtained from last-generation realistic potentials, and used a realistic model for the nuclear weak current, 
including one- and two-body terms. The vector part of this current was shown to provide an excellent description of 
the np radiative capture cross section for neutron energies up to fOO MeV [11], while the axial part was constrained 
to reproduce the Gamow- Teller matrix element in tritium /3-decay [12]. The Nakamura et al. studies have played 
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an important role in the analysis and interpretation of the Sudbury Neutrino Observatory (SNO) experiments [13], 
which have established solar neutrino oscillations and the validity of the standard model for the generation of energy 
and neutrinos in the sun [14]. 

In the present work, we use the same theoretical framework as the authors of Refs. [11, 12], but include refinements 
in the modeling of the weak current — which however, as shown in Sec. V, will turn out to have a minor impact on 
the predicted cross sections — and extend the range of neutrino energies up to 1 GeV. While the theoretical approach 
is essentially the same, nevertheless the way in which the calculations are carried out in practice is rather different 
from that used in those earlier papers, which relied on a multipole expansion of the weak transition operators, and 
evaluated the cross section by summing over a relatively large number of final two-nucleon channels states. In contrast, 
we evaluate, by direct numerical integrations, the matrix elements of the weak current between the deuteron and the 
two-nucleon scattering states labeled by the relative momentum p (and in given pair spin and isospin channels), 
thus avoiding cumbersome multipole expansions. Differential cross sections are then obtained by integrating over 
p (and summing over the discrete quantum numbers) appropriate combinations of these matrix elements, i.e. by 
calculating the weak response functions. The techniques developed here for the deuteron should prove valuable when 
we will attempt the Green's function Monte Carlo calculation of these response functions (or rather, their Laplace 
transforms [15]) in A > 2 nuclei. 

This paper is organized as follows. In Sec. II and Appendix A we present the neutrino and antineutrino differential 
cross sections expressed in terms of response functions, while in Sec. Ill we provide a succinct description of the neutral 
and charge-changing weak-current model. In Sec. IV we outline the methods used to obtain the two-nucleon bound 
and continuum states, and discuss the numerical evaluation of the response functions. A variety of results for the 
neutral current processes 2 H{vi,vi)pn and 2 H(l^, v{)pn, and charge- changing processes 2 H(j/ e , e~)pp and 2 H(l7 e , e+)nn 
are presented in Sec. V, including the sensitivity of the calculated cross sections to (i) interaction effects in the final 
states, (ii) different short-range behaviors of the two-body axial weak currents, and (iii) different potential models 
to describe the two-nucleon bound and continuum states. In order to illustrate the effects of nuclear structure, we 
compare these cross sections to those obtained in a naive model in which the deuteron is taken to consist of a free 
proton and neutron (the free nucleon cross sections are listed for reference in Appendix B). Concluding remarks and 
an outlook are given in Sec. VI. 



II. INCLUSIVE NEUTRINO SCATTERING OFF DEUTERON 



The differential cross section for neutrino {v) and antineutrino {v) inclusive scattering off deuteron, specifically the 
processes 



n + d — > vi+p + n , vi+d — > vi+p + n 
induced by neutral weak currents (NC), and the processes 

ui+ d — > l~ + p + p , T>i+ d — > l + + n + n 
induced by charge-changing weak currents (CC), can be expressed as 
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where G—Gf for the NC processes and G=Gf cos 6c for the CC processes, and the — (+) sign in the last term is 
relative to the v (v) initiated reactions. Following Rcf. [12], we adopt the value Gf = 1.1803 x 10~ 5 GcV~ 2 as obtained 
from an analysis of super-allowed + — > + /3-decays [16] — this value includes radiative corrections — while cos 9c is 
taken as 0.97425 from [17]. The initial neutrino four-momentum is k 11 — (e, k), the final lepton four momentum is 
k^' = (e', k'), and the lepton scattering angle is denoted by 9. We have also defined the lepton energy and momentum 



transfers as 



The Fermi function F(Z, k') with Z 
the charge-raising reaction, 



and q = k — k', respectively, and the squared four- momentum transfer as Q 



oj 2 >0. 



2 accounts for the Coulomb distortion of the final lepton wave function in the 



F(Z, k 1 ) = 2 (1 + 7 ) (2 k 1 r d f ^ 2 exp (tt y) 



r(i + 2 7 ) 
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it is set to one otherwise. Here y = Z ae'/k', T(z) is the gamma function, r<j is the deuteron radius (r^ = 1.967 fm), 
and a is the fine structure constant. Radiative corrections for the CC and NC processes due to bremsstrahlung and 
virtual photon- and Z-exchangcs have been evaluated by the authors of Rcfs. [18, 19] at the low energies (<~ 10 MeV) 
relevant for the SNO experiment, which measured the neutrino flux from the 8 B decay in the sun. These corrections 
are neglected in the present work, since its focus is on scattering of neutrinos with energies larger than 100 MeV. 
We are not (or not yet, at least) concerned with providing cross section calculations with % accuracy in this regime. 
Lastly, the nuclear response functions are defined as 
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where | d, M) and |/) represent, respectively, the initial deuteron state in spin projection M and the final two-nucleon 
state of energy Ef, and is the deuteron rest mass. The three-momentum transfer q is taken along the £-axis (i.e., 
the spin-quantization axis), and j^(q, w) is the time component (for /j = 0) or space component (for \i = x,y,z) of 
the NC or CC. 

The expression above for the CC cross section is valid in the limit e' ~ k', in which the lepton rest mass is 
neglected. At small incident neutrino energy, this approximation is not correct. Inclusion of the lepton rest mass 
leads to changes in the kinematical factors multiplying the various response functions. The resulting cross section is 
given in Appendix A. 



III. NEUTRAL AND CHARGE-CHANGING WEAK CURRENTS 

We denote the neutral and charge- changing weak currents as and Jq C , respectively. The former is given by 

f NC = -2 sm 2 9 w 3^ s + (1 - 2 sin 2 ^) j* z + jf , (3.1) 

where 9 W is the Weinberg angle (sin 2 ^ = 0.2312 [17]), and denote, respectively, the isoscalar and isovector 
pieces of the electromagnetic current, and denotes the isovector piece of the axial current (the z on the isovector 
terms indicates that they transform as the ^-component of an isovector under rotations in isospin space). Isoscalar 
contributions to associated with strange quarks are ignored, since experiments at Bates [20] and JLab [21] have 
found them to be very small. 

The charge-changing weak current is written as the sum of polar- and axial-vector components 

f cc = jt+j£ , 3±=jx±iJv (3-2) 

The conserved-vector-current (CVC) constraint relates the polar-vector components j£ of the charge-changing weak 
current to the isovector component j0j z of the electromagnetic current via 

[T a , j^ z ] =ie azb j£ , (3.3) 

where T a are isospin operators. We now turn to a discussion of the one- and two-body contributions to the NC and 
CC. 
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A. One-body terms 

The isoscalar components of the one-body electromagnetic current are given by 

Gf(Q 2 ) n .2Gf f (Q 2 )-G|(Q 2 ) 
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and the corresponding isovector components of z are obtained by the replacements 

G|(Q 2 ) — ► G v E {Q 2 ) Tl , z , G S M (Q 2 ) — ► G^(g 2 )r JjZ , 



(3.4) 

(3.5) 
(3.6) 
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where G^/ y and Gf[ V are the isoscalar/isovector combinations of the proton and neutron electric (E) and magnetic 
(M) form factors, and pi are the position and momentum operators of nucleon i, cti and t, i2 are its Pauli spin 
and isospin matrices, and m is the nucleon mass (0.9389 GeV). Note that we have decomposed j 7i g and j 7iZ into 
transverse (_L) and longitudinal (||) components to the momentum transfer q, and have used current conservation to 
relate the latter to the isoscalar and isovector charge operators j® s and j® z . The isovector components of the axial 
weak neutral current j^ 5 are given by 



G A (Q 2 ) 
4 m 

G A (Q 2 ) 



Ti, <Ji 



[Pi : 



(3.8) 



er, ; c 



4 m 2 



o-i [p- , c lq ' r *] + - [(ffi-pOpi, e l 



■ q [Pi , e tq ' r *] + - -q [<Ti • Pi 



3 J qr; 



■i q x pi e 



(3.9) 



where G^ is the nucleon axial form factor, and [. . . , . . .], denotes the anticommutator. The operators above include 
terms of order (v/c) 2 in the non-relativistic expansion of the single-nucleon covariant currents. These have been 
neglected in the study of Rcf. [12]. The proton and neutron electromagnetic and nucleon axial form factors are 
parametrized as 



G p E {Q 2 ) 

G p M (Q 2 ) 
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G D (Q 2 ) , G n E {Q 2 ) 



Q 2 G D (Q 2 ) 



4m 2 1 + Q 2 /™ 2 ' 
H P G D (Q 2 ) , G n M (Q 2 )=fi n G D (Q 2 ) , 
1 ~ , n2 N _ 9a 



(1 + Q 2 /A 2 ) 2 



G A {Q 2 ) = 
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(3.10) 
(3.11) 
(3.12) 



from which the isoscalar and isovector combinations are obtained as G EM = G P E M ± G E M . The proton and neutron 
magnetic moments are [i v — 2.793 and /j, n — —1.913 in units of nuclear magnetons (n.m.), and the nucleon axial-vector 
coupling constant is taken to be g A — 1.2694 [17]. The values for the cutoff masses A and used in this work are 
0.833 GeV and 1 GeV, respectively. The former is from fits to elastic electron scattering data off the proton and 
deuteron [22], while the latter is from an analysis of pion electroproduction [5] and neutrino scattering [6-8] data. 
Uncertainties in the Q 2 dependence of the axial form factor, in particular the value of A^, could significantly impact 
predictions for the neutrino cross sections under consideration. As mentioned earlier, recent analyses of neutrino 
quasi-elastic scattering data on nuclear targets [4] quote considerably larger values for Aa, in the range (1.20-1.35) 
GeV. 

The polar- vector (j±) and axial- vector (j± 5 ) components of the charge- changing weak current are obtained, respec- 
tively, from jUj z and by replacing 



1~i,z/2 > Ti t ± = (Ti >x ± Ti. y )/2 . 



(3.13) 
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However, in the case oi j± 5 , in addition to the terms entering Eqs. (3.8)-(3.9), we also retain the induced pseudoscalar 
contribution, given by 



PS) = - G PsiQ2) n. ± q » a, ■ q e*™ , 
2 m m M 

where the induced pseudoscalar form factor Gps is parametrized as 



(3.14) 



(3.15) 



This form factor is not well known [23]. The parameterization above is consistent with values extracted [24, 25] from 
precise measurements of muon-capture rates on hydrogen [26] and 3 He [27], as well as with the most recent theoretical 
predictions based on chiral perturbation theory [28]. This contribution vanishes in NC-induced neutrino reactions. 

B. Two-body terms 

Two-body terms in the neutral and charge-changing weak currents have been discussed in considerable detail in 
Refs. [29-31] (and references therein). We list the terms included in the present study — i.e., the subset of those derived 
in the above references expected to give the dominant two-body contributions to the processes of interest here — in 
the following two sub-sections for clarity of presentation and future reference in Sec. V. Unless stated otherwise, they 
are given in momentum space, and configuration-space expressions follow from 

0(q)= / III ^^(k.+k.-qK^+^/V*--^-^ 

J\n JKi Jk, JK, 



where kj = p[ — pi and Kj = (p^ + Pi)/2, pi and Pj are the initial and final momenta of nuclcon i, and 

dp 



P , (2^) 3 " 

These configuration-space operators are used in the calculations reported below. 



(3.17) 



1. Two-body vector terms 

The two-body isovector current operator j 7 ,z(y) consists of pseudoscalar- and vector-meson (referred to as 7r-likc 
and p-like) exchange, and A-excitation terms, 



C—7T, p, A 



(3.18) 



The 7r-like and p-like exchange currents read: 



hAV'^) = iGsiQ )(T; X Tj) z V n (k 
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h,z(ij;p) = -iG v E {Q 2 ){n x Tj) 2 



+(<Ti x kj) a-j ■ (kj x kj) + (<Xj x kj) <Tj • (kj x kj) - v pS (kj 



(kj - kjXo-j x kj) ■ (a-j x kj) 



^ ^3 

kf~k 2 



(3.20) 



where 



v n (k) = v° T {k) - 2 v tT (k) , Vp (k) = v aT (k) + v tT {k) , v p s(k) = v T (k) 



(3.21) 
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and 

/•OO 

v T (k) = 4tt / r 2 dr j (kr)v T (r) , (3.22) 
Jo 

A /"OO 

yCTT ( fc ) = F/o ^ bo(fcr)-l]^ r (r) , (3.23) 

v tT {k) = -j~2 I r 2 drj 2 (kr)v tT (r) . (3.24) 
* Jo 

Here v T (r), u°" r (r), t> tT (r) are the isospin-dependent central, spin-spin, and tensor components of the two-nucleon 
interaction (the AV18 in the present study), and ji(kr) are spherical Bessel functions. The factor jo(kr) — 1 in the 
expression for v aT (k) ensures that its volume integral vanishes. In a one-boson-exchange (OBE) model, in which the 
isospin-dependent central, spin-spin, and tensor interactions are due to ir- and p-meson exchange, the functions v v (k), 
v p (k), and v p s(k) simply reduce to 

~+ -^kwvk ' (3 - 25) 

Vp{k) ^ 4m2 k^+my {d - 2b) 

o h 2 Jk) 

-► - (3-27) 

where and m p are the meson masses, f n , and g p and k p are the pseudovector irNN, and vector and tensor pNN 
coupling constants, and the hadronic form factors are parameterized as 

h °W = l^TW> a = Tr,p. (3.28) 

While the AV18 interaction is not a OBE model, nevertheless the effective propagators v v (k), v p (k), and v p s(k) 
projected out of its v T (k), v (JT (k), and v tT (k) components are quite similar to those listed above with cutoff masses 
in the range (1.0-1.5) GeV. We note that the 7r-like and p-like currents with the v n (k), v p (k), and v p s(k) defined in 
Eq. (3.21) satisfy by construction the current conservation relation with the AV18 t, or, and tr interaction components 
(for a discussion of the issue of current conservation in relation to the momentum-dependent terms of the AV18, see 
Ref. [31]). 

The isovector A-excitation current is written in configuration space as (for a derivation based on a perturbative 
treatment of A-isobar degrees of freedom in nuclear wave functions, see Ref. [29] ) 

• Hz {ij\ A) = -i 2 ^m iQ m A ) 6iq ' r ' Kiv(y) qxS, T i>z + adjoint] , (3.29) 

where S and T are spin- and isospin-transition operators converting a nucleon into a A isobar and satisfying the 
identity 

S f • A S • B = \ A • B - % -rr • (A x B) , (3.30) 

v AN(ij) is the NN to AN transition potential, 

VAN(ij) = [ vl r N ( rij ) S, • aj + vX N { rij ) S% N } T, • Tj , (3.31) 

Sfj N is the tensor operator obtained by replacing <Xj with Sj, the regularized spin-spin and tensor radial functions 
w A7v( r ) an d "Ajv( r ) are defined as 

< fc,< p) _M!*( 1+ | + »)£( 1 _.-^)", (3.33) 
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and x — rn^r, f* is the ttA^A coupling constant (/* = 2.19 from the width of the A), and the parameter in the 
short-range cutoff function is taken as A = 4.29 (from the AV18). Finally, the 7./VA electromagnetic transition form 
factor G 7 na is parameterized as 

G JNA (Q 2 ) = - M ^ = , (3.34) 



(l + QVAi^yi + QVA^ ' 



where the transition magnetic moment /i^na is 3 n.m., as obtained in an analysis of 7 TV data in the A-resonance 
region [32]. This analysis also gives Aa,i=0.84 GeV and Aa,2=1-2 GeV. 

The two-body isoscalar current operator j 7l s(«i) considered in the present study only includes the contribution 
associated with the pn-f transition mechanism, 

h,s{ij) = h,s(ij; pn)+i^ j , (3.35) 

where 

f g h (k ) h (k ) 

h,s(ij;pn) = - iG p^(Q^9ptt~/^^ Ti-rj k / + l m2 fc / + 3 m2 (ki x k io- -k, , (3.36) 

The combination of coupling constants gp^f^gp is taken as 1.37, and the cutoff masses A^ and A p as 0.75 GeV 
and 1.25 GeV, respectively, from a study of the deuteron magnetic form factor [33]. The Q 2 dependence of the 
electromagnetic transition form factor G p7r7 (Q 2 ) is modeled by using vector- meson dominance, 

G P ^(Q 2 ) = 1 - J 2/m 2 , (3.37) 

where m u is the w-meson mass. 

The two-body isovector and isoscalar electromagnetic charge operators j® z and j® s consist of terms associated 
with 7r-like and p-like exchanges 

3*,z/s(ij)= E [j°, z /s(ir,c)+i^j] , (3.38) 



C—TT, p 



where 



and 



F v (O 2 ) 

jjAV' 7r ) = 1 2m v ^ k o) ( a i ■ q) (°j ' k i) ' ( 3 - 39 ) 

3°, x (ir, P) = Mkj) x q) • (trj xkj), (3.40) 

Zi Til 



fi,s(ir, = F ^P- Ti ■ Tj V v (kj) (<Tj • q) (*j ■ k,) , (3.41) 



Ff(Q 2 ) 



with v n {k) and v p {k) as defined in Eqs. (3.21). The nucleon electromagnetic Dirac and Pauli form factors F^ V and 

s/v 

F 2 are obtained from 

F siv m = G S J V (Q 2 )+^G% v {Q 2 ) (3 43) 

F s /V{Q2) = gT{q 2 )^g s J v (Q 2 ) (3 44) 

with 17 = Q 2 /(4m 2 ). 

The polar-vector components j± of the charge-changing weak current j£ c are obtained from z via CVC, which 
implies the replacements 

{Ti X Tj) z > {Ti X Tj)± = {Ti X Tj) x ±i{TiX Tj) y (3.45) 
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in Eqs. (3.19)-(3.20), 

T t . z /2 — > T i>± = (T i;X ± i Ti >y )/2 (3.46) 

in Eq. (3.29), and the replacement (3.13) in Eqs. (3.39)-(3.40). Only the transverse components (perpendicular to q) 
of the vector part of the NC and CC are explicitly included in the calculations to follow. Their longitudinal components 
have already been effectively accounted for by the replacement in Eq. (3.6) (and the similar one for the isovector terms). 
Lastly, we note that in the study of Rcf. [12] the p-meson exchange and pir transition contributions to the two-body 
vector current, and tt- and p-exchange contributions to the two-body vector charge have been neglected. Furthermore, 
the 7r-exchange and A excitation currents are regularized by introducing a monopole form factor (A ff = 4.8 fm -1 ), 
which naturally leads to a different short-range behavior of these currents than obtained here. 



2. Two-body axial terms 



TABLE I: Contributions to the GT matrix element in tritium /3-decay. The one-body (1-b) NR and RC contributions are, 
respectively, from the leading and 1/m 2 terms in Eq. (3.9); the two-body (2-b) contributions are from Eqs. (3.48)-(3.51). Set 
I (II) corresponds to the cutoff choices A^ = A p = 1.2 GeV (A-n- = 1.72 GeV and A p = 1.31 GeV). The N to A axial coupling 
constant g* A for each set is obtained by fitting the experimental value of the GT matrix element, given by 0.955 ± 0.002 [25]. 





Set I 


Set II 


1-b (NR) 


+0.9213 


+0.9213 


1-b (RC) 


-0.0085 


-0.0085 


2-b (tt) 


+0.0078 


+0.0123 


2-b (p) 


-0.0042 


-0.0055 


2-b (pn) 


+0.0123 


+0.0196 


2-b (A) 


+0.0263 


+0.0159 


9*a/9a 


0.614 


0.371 



The axial parts of the neutral and charge-changing weak current operators consist of contributions associated with 
tt- and p-meson exchanges, the axial p7r transition mechanism, and a A excitation term 



C—7T. p, p7T. A 



(3.47) 



where the isospin component a is either z for NC or ± for CC. The ir- and p-meson exchange and p7r transition axial 
currents read, respectively, 



G A (Q 2 ) ft hl{k 3 ) 



2 m ml k 2 + ml 



(n x Tj) z a-i x kj - Tj >z (q + 2icr i x Kj) 



<Tj ■ k, , 



(3.48) 



Ga(Q 2 ) 9 2 P (l + k p ) 2 h^kj) 



4 m 2 k 2 j + to 2 , 



2to 

\°i x ("j x k j)] x kj- 



+ L 



3,z 



(n x Tj) z q <Ji ■ (a-j x kj) + 2i (crj xk 3 )x 
(a-j x kj) x kj ■ — 2i [o-j x (a-j x kj)] x Kj 



G A (Q 2 ) 2 h p (h) K(kj) 



TO 
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P k 2 + TO 2 k 2 + TO 2 



(TiXT 3 ) z (l + K p )cr i xk i -2iK l 



CTj • kj , 



while the A excitation axial current is obtained from [29] 



j^;A) 



2 (m — toa) 

where the (unknown) A^ to A axial form factor is parameterized as 

G* A (Q 2 ) = ^ , • 

(1 + Q 2 /A 2 4 ) 2 



^Aiv(^) S i T i,z + adjoint 



(3.49) 



(3.50) 



(3.51) 



(3.52) 
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The charge-changing axial currents follow by replacing the isospin operators as in Eqs. (3.13) and (3.45)-(3.46). 
The values for the it- and p-meson coupling constants are taken from the CD-Bonn one-boson-exchangc model [34], 
/2/(4tt) = 0.075, 3p/(47r) = 0.84, and n p = 6.1, while two different sets of cutoff masses A v and A„ are used in the 
present work: A x = A p = 1.2 GeV (Set I) in line with the cutoff masses extracted from the 7r-like and p-like exchanges 
associated with the AV18 model; A,, = 1.72 GeV and A p = 1.31 GeV (Set II) from the CD-Bonn model. In the N to 
A axial current, the Q 2 dependence of the form factor is taken to be the same as that of the nucleon; however, the 
value for the transition axial coupling constant g* A is determined by fitting the Gamow- Teller (GT) matrix element of 
tritium /3-decay [25] in a calculation based on trinucleon wave functions corresponding to the AV18/UIX Hamiltonian 
and the present model for the axial current. The values corresponding to Set I and II of cutoff masses are listed in 
Table I. 

Finally, in the present study the axial charge operator is taken to include only the pion-exchange term, whose 
structure and strength are determined by soft-pion theorem and current algebra arguments [35] 

3a 5 (i3)=3a 5 (ir,K) + i^3 , (3-53) 

where 

i- M (y ! *) - "< fa >< rj)a *i ■ k, , (3.54) 

is pion decay amplitude (F v = 93 MeV), and the Q 2 dependence of the associated form factor is assumed to be 
the same as in the nucleon. We conclude by noting that the model described above for the two-body axial charge and 
current operators is essentially the same of that used in Ref. [12], apart from differences in the values of the cutoff 
masses for the hadronic form factors of the meson exchange terms, and a different treatment of the A excitation 
current. However, it is important to stress that both here and in Ref. [12] the two-body axial currents are constrained 
to reproduce the experimental tritium /3-decay rate. 



IV. CALCULATION 



The two-body scattering- and bound-state problems are solved in momentum space with the methods discussed in 
Ref. [36], which facilitates calculations with a non-local potential such as CD Bonn. We briefly summarize them in 
the next two sub-sections for clarity. In the last sub-section we discuss the calculation of the weak current matrix 
elements, response functions, and cross sections. 



A. The scattering-state problem in momentum space 



In the case of scattering (setting aside the treatment of the Coulomb interaction for the time being), we solve for 
the if-matrix in channel J ST (hereafter, L is the relative orbital angular momentum, S and T are the total spin and 
isospin, and J is the total angular momentum, and (— i) L + s + T = — l) 



V JST 
V L'L 



(p';p) + 



4p 



dkk 



2 E< 

L" 



V 



Kjf, T L (k;p) 



(4.1) 



where p is the two- nucleon reduced mass, V denotes a principal- value integration, and v J L f^{p' : p) are the p-space 
matrix elements of the potential, projected in channel J ST [36]. We should note the presence of the somewhat 
unconventional phase factor i L ~ L included in the matrix elements vjf£ (p';p) [36], which makes the states used here 
differ by a factor i L from those usually adopted in nucleon-nucleon scattering analyses. The integral equations (4.1) are 
discretized, and the resulting systems of linear equations are solved by direct numerical inversion. The principal- value 
integration is eliminated by a standard subtraction technique [37]. Phase shifts in channel J ST are easily obtained 
from the on-shell S'-matrix related to the (on-shell) if-matrix by 



S JST (p)= [l + 2ip P K JST (p;p)] 1 [1 -2ippK JST {p;p) 
while r-space wave functions follow from 



*#2(r;j>) 



j(pr) + — r dkk 2 j{kr)^—K JST {k-p) 

TT JO p K 



J L'L" L 



l + 2ippK Jbl (p;p) 



(4.2) 



(4.3) 



L"L 
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where the matrix [j(pr)}L>L = &VL 3l{p t ) has been introduced for convenience. The (complex) radial wave functions 
Zjflir) behave in the asymptotic region r — > oo as 



z J L %\r;p) * - [6 LIL hf{pr) + h^(pr)S^ (p) 



(4.4) 



where the functions h L ' (pr) are defined in terms of the regular and irregular (n^) spherical Bessel functions as 

h^ 2 \y)=j L {y)±in L (y) . (4.5) 

In the calculation of the response functions that follows, scattering wave functions with incoming-wave boundary 
condition (— ) are required. These are written as 



(4.6) 



where Xm s anc ^ ^Mt are two-nucleon spin and isospin states, respectively, y^fsj are standard spin-angle functions, 



Z J lsm s (p) = ^2( lm l,SM s | JMj) 



(4.7) 



and (LMl, SMs \ JMj) are Clebsch-Gordan coefficients. Note that the wave function in Eq. (4.6) retains interaction 
effects only in channels with J < J max , and reduces to plane waves for J > J max . 

When the Coulomb interaction is present, we use the method developed originally in Ref. [38], which allows us to 
solve the pp scattering problem in momentum space [39] . It consists essentially in separating the potential into short- 
and long-range parts v$ and vl 1 where vl only includes the Coulomb potential vc and vs includes, in addition to 
vc, the nuclear potential v. Then the standard momentum-space technique outlined earlier can be used to solve the 
problem with v$, and the corresponding radial wave functions behave as 



(4.8) 



where Sg S1 is the S'-matrix in this case (with T = 1), and the are normalization constants. The wave functions 
z s S L'L should match smoothly those relative to vs + vl, which behave asymptotically as 



z J L ?£{r-p) ~ \ \5 L , L h { l\i,pr) + h^^pr) S J L ?l(p) 



(4.9) 



where 



h ( L 2) (£, y) = [F L & y) T G L (^ y)]/y, t = a »/p , 



(4.10) 



and Fl and Gl are the regular and irregular Coulomb functions. Carrying out the matching for the functions 
and their first derivatives leads to a relation between Sg S1 and S JS1 and a determination of the normalization 
constants [39]. Finally, pp scattering wave functions with incoming- wave boundary conditions are written as in 
Eq. (4.6) with T, Mr = 1,1 and the replacement 



y JS\* 

Z L'L 



(r;p) 



- Z L'L 



{r;p) , 



(4.11) 



where ol is the Coulomb phase shift, 



CT£=arg[T(L + l + iO] • 



(4.12) 



Hence, Coulomb interaction effects are retained only in channels with J < J ma x, and are ignored in those with 
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B. The bound-state problem in momentum space 

The deuteron wave function is written in r-space as 

fe(r)= ]T l L u L (r)yl[ 1 {v)r l l, (4.13) 

L=0,2 

and the radial wave functions Ul{t) (L = 0, 2) follow from 

2 f°° 

u L {r) = -j dpp 2 j L (pr)u L (p) . (4.14) 
f Jo 

The p-space wave functions Ul{p) are obtained from solution of the homogeneous integral equations 

u L ( P )=r 277^-/ dfcfc2 E v™{p-k)u L ,(k) . (4.15) 

Here, is the nuclear potential in the JST = 110 channel, and E d is the deuteron energy (E d = —2.225 MeV). 
We again note the unconventional phase i L in Eq. (4.13). 

C. Matrix elements, response functions, and cross sections 

The deuteron wave function in Eq. (4.13) is written, for each spatial configuration r, as a vector in the spin-isospin 
space of the two nucleons, 

^ M (r) = ^< l) (r)|n), (4.16) 

n=l 

where | n) = (p t)i (n |) 2 , (n t)i (p t)2, • • • , ( n i)i (p D2 an d ipffl are the components of m this basis. In the 
NC-induced processes, the scattering wave function in Eq. (4.6) is expanded on the same basis; however, in the CC- 
induced processes the pp or nn scattering wave functions are expanded on the (spin only) basis | m) , , and 
11 for pp or nn. Matrix elements of the weak current operators are written schematically as 

(f\0\d,M)=J dr^V^'WOn'.nW^W . (4-17) 

where the momentum- and energy-transfer dependence is understood. The spin-isospin algebra is performed exactly 
with techniques similar to those developed in Ref. [42], while the r-space integrations are carried out efficiently by 
Gaussian quadratures. Note that no multipolc expansion of the transition operators is required. When momentum 
operators are present, they are taken to act on the right (deuteron) wave function. For example, the one-body axial 
charge operator is written as 

O(r) ^M(r) — ► [e iqT/2 <Ti • (-2iV + q) r M + 1^2] 1> M (r) , (4.18) 

and the derivatives are evaluated numerically as 

where e a is a unit vector in the a-direction, and S is a small increment. Once the matrix elements have been computed, 
response functions are evaluated (in the lab frame) via 



= *E E I i^\^ + rn d -E + -E_)f s ab M ^ M {w), 



(4.20) 



M SM S ,T' 

with 

/ a 5 fc Ms '™ T;M (q,p) = (c L ,p;SM s ,TM T \O a (c l ,uj)\d,M)(c l ,p;SMs,TM T \O b (ci;u;)\d,Mr , (4.21) 
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where |q, p; SMg, TMt) represents the final two-nucleon scattering state with total momentum q and relative mo- 
mentum p, raj is the deuteron rest mass, and E± are the nucleons' energies in the final state, 



E± = v/(q/2 ± p) 2 + m 2 . 



(4.22) 



The factor 1/2 in Eq. (4.20) is to avoid double-counting the contribution of the final states (which are anti- 
symmetrized), and the pair isospin T assumes the values T = 0, 1 with M T = for NC processes, and T = 1 
with Mt = 1 or —1 for CC processes. The 5-function is integrated out, and 



i^^/! dxp 



M SM S ,T 



p + xq/2 ^ p — x q/2 



E 4 



f SMs,TM T ;M [q ^ x) 



(4.23) 



where x = q • p, and the magnitude p of the relative momentum is fixed by energy conservation. This magnitude 
depends on q, lo, and x. However, in order to reduce the computational effort, the scattering states entering the 
product of matrix elements f a b are obtained at the energy 



4(p 2 +m 2 ) = (uj + m d ) 2 -q 2 



(4.24) 



which only depends on q and uj. Lastly, Gauss points (~ 50) are used to perform the x-integration accurately. 
Extensive and independent tests of the computer programs have been completed successfully. 

Total cross sections are obtained by direct integration of Eq. (2.3) by evaluating the differential cross sections on a 
grid of Gauss points in e' (the lepton final energy) and 9 (its scattering angle). There are kinematical constraints on 
the allowed values for e' and 0, which follow from the requirement p 2 > 0: 



ATo 2 o-^f , \ { i —\ - m d (e- 6th) +m t {mi+ 2m) 
e J e' 2 - mf cos 9 > (e + m d ) (e — e) , e = 



where e t h is the threshold energy for the initial neutrino (e > e t h), 

(mi + 2m) 2 - m 2 d 



eth 



2m d 



(4.25) 



(4.26) 



mi is the rest mass of the final lepton (mi — in the NC case), and m = (m p + m n )/2 for NC reactions or m = m p 
(m n ) for charge-raising (charge- lowering) reactions. These kinematical constraints imply: 



mi<e' < e'_ for - 1 < cos 9 < , 
mi < e' < e' + for < cos 9 < 1 , 



where the limits e'± are defined as 



e±^e 2 - (l-f3 2 cos 2 9) (e 2 +m 2 ^ 2 cos 2 9) 
1 - B 2 cos 2 9 



1 + m d /e 



In the case of NC reactions (m; = 0), they are simply given by 



< e' < 



(3 cos 9 



for 



1< cos 9 < 1 . 



(4.27) 
(4.28) 



(4.29) 



(4.30) 



V. RESULTS 



Cross section values obtained with the AV18 interaction and the one- and two-body terms in the electroweak 
current discussed in Sec. Ill are listed in Tables II-IV for initial neutrino energies in the range (5-1000) MeV. The 
two-body axial currents are those corresponding to Set I in Table I. The two-nucleon (NN) scattering states are 
written as in Eq. (4.6): they include interaction effects in channels with J < J max = 5 and reduce to spherical Bessel 
functions (i.e., plane waves) for J > J max . The relative kinetic energy T = 2 (p 2 + m 2 ) 1 / 2 — 2 m, where p is defined 
in Eq. (4.24), changes over a wide range as the initial neutrino energy increases up to 1 GeV and the final lepton 
energy and scattering angle vary over the allowed kinematical regions: at e = 50 MeV, < T < 48 MeV; at e = 500 
MeV, < T < 445 MeV; and at e = 1000 MeV, < T < 819 MeV. These relative energies (at the larger values of 
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TABLE II: Total cross sections in cm 2 for the NC- and CC-induced processes on the deuteron as function of the initial neutrino 
energy e, obtained with the AV18 potential and the inclusion of one- and two-body terms in the weak current. The number in 
parentheses, "—a;", denotes 10~ x ; for instance an entry like 9.561(-44) stands for 9.561xlCP 44 cm 2 . 



(MeV) 



5 
10 
20 
30 
40 
50 
60 
70 
80 
90 
100 



9. 561 (-44) 
1.104(-42) 
6.965(-42) 
1.833(-41) 
3.555(-41) 
5.892(-41) 
8.839(-41) 
1.240(-40) 
1.657(-40) 
2.131(-40) 
2.657(-40) 



-NC 



9.363(-44) 
1.053(-42) 
6.285(-42) 
1.568(-41) 
2.885(-41) 
4.546(-41) 
6.495(-41) 
8.699(-41) 
l.lll(-40) 
1.369(-40) 
1.640 (-40) 



-CC 



3.427(-43) 
2.680(-42) 
1.547(-41) 
4.058(-41) 
7.995(-41) 
1.348(-40) 
2. 338 (-40) 
2. 949 (-40) 
4.036(-40) 
5.320(-40) 
6.631(-40) 



-CC 



2. 831 (-44) 
1.242(-42) 
9.562(-42) 
2.508(-41) 
4.685(-41) 
7.403(-41) 
1.057(-40) 
1.409(-40) 
1.790 (-40) 
2.191(-40) 
2. 606 (-40) 



e) are well beyond the range of applicability of all modern realistic interactions, which are typically constrained to fit 
NN scattering data up to pion production threshold (T ~ 150 MeV). This is also the case for the AV18 of course, 
although it is known [40] that this interaction reproduces quite well phase shifts (at least in those channels where 
inelasticities are small) up to T < 300 MeV. 

As an additional caveat, we note that the present theory cannot describe the inclusive cross section in the pion- 
production region, for example the A-excitation peak region, since no mechanisms for (real) single- or multi-pion 
production are included in it. However, it does reproduce quite well the observed d(e, e') inclusive cross section in 
the quasi-elastic peak region at intermediate values of the three-momentum transfer. This is illustrated in Fig. 1 
(a-d), where the longitudinal and transverse response functions Rl and Rt obtained at Bates [41] by Rosenbluth 
separation of (e, e') data at momentum transfers of 300 MeV and 500 MeV are compared with theory. In these 
figures, we show separately the response functions calculated with an electromagnetic current including one-body 
only and (one-l-two)-body terms, as well as those obtained by replacing the fully interacting NN states of Eq. (4.6) 
with plane waves (curves labeled by PW). Two-body terms in Rl give negligible contributions, those in Rt lead to 
an increase of the transverse strength over the whole quasi-elastic region, which amounts to a few % at the top of 
the peak, but becomes sizable (relative to the one-body response) as the energy transfer u increases well beyond the 
quasi-elastic peak energy w qc = (q 2 + to 2 ) 1 / 2 — m. At these momentum transfers, the quasi-elastic and A peaks in 
Rt, the latter at cja = (q 2 + m 2 ^) 1 / 2 — m (toa = 1232 MeV), are well separated — note, however, the rise seen in 
the data at q = 500 MeV and the highest w's, presumably due to (transverse) strength creeping in from the A-peak 
region. Interaction effects in the NN continuum states are important, particularly at low momentum transfers and/or 
for energy transfers close to the threshold for deuteron breakup. However, at the larger q-values plane-wave states 
provide response functions in the quasi-elastic region, which are fairly close to those predicted by the exact scattering 
states. Finally, we note that at q = 500 MeV and quasi-elastic energies theory over-predicts the measured longitudinal 
response. As a consequence, the total integrated longitudinal strength — the Coulomb sum rule — obtained from these 
data [42] is smaller than calculated. On the other hand, there is excellent agreement between the theoretical and 
measured Coulomb sum rules at q = 300 (and 400) MeV [42]. 

In Table II the columns labeled t/;-NC and I7;-NC refer to the NC-induced reactions in Eq. (2.1), those labeled 
v e -CC and V e -CC refer to CC-induced reactions in Eq. (2.2), and the initial neutrino energy is between 5 MeV (close 
to threshold) and 100 MeV. In this energy range the cross sections change rapidly, by 3-4 orders of magnitude, and 
interaction effects in the final scattering states are important. Two-body terms in the vector and axial pieces of the 
weak current increase the one-body cross section typically by a 2-3 % for both the NC- and CC-induced reactions, in 
agreement with the results of Rcf. [12]. 

There are differences between the present calculations and those of Nakamura et al. [12] — mostly having to do with 
the model for the weak current — which, however, only lead to small numerical differences in the predicted cross section 
values, as shown below. As we have already remarked in Sec. IIIB, the authors of that work ignore the relativistic 
corrections proportional 1/m 2 in the one-body axial current (3.9), and use a different short-range parametrization for 
the two-body vector and axial currents than adopted here. In addition, the cutoff masses entering the nucleon form 
factors have slightly different values from those listed in Sec. Ill A. In order to have a more meaningful comparison 
with the results of that work, we have carried out a calculation of the NC- and CC-induced cross sections at three 
representative initial neutrino energies, in which we have removed the relativistic correction in the one-body axial 
current and have changed the cutoff mass values in the nucleon form factors so as to match those used in Ref. [12]. 
Inspection of Table V shows that the two calculations are typically within less than 1% of each other. This level of 
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FIG. 1: (color online) The electromagnetic responses: longitudinal at q — 300 MeV (a); longitudinal at q = 500 MeV (b): 
transverse at q = 300 MeV (c); transverse at q = 500 MeV (d). These are obtained with the AV18 potential and the inclusion of 
one-body (dashed line) and (one+two)-body (solid line) terms in the nuclear electromagnetic charge operator, and are compared 
to data. Also shown are the results obtained with plane-wave (PW) final states. 



agreement should be viewed as satisfactory, given the different ways in which the two calculations are carried out 
in practice. The authors of Ref. [12] rely on a multipolc expansion of the cross section, whereas we compute the 
matrix elements entering the various response functions by direct numerical integrations, which avoid the need for 
introducing (cumbersome) multipole expansions of the weak transition operators. Of course, the present calculations 
are computationally intensive: evaluation of the NC cross sections requires about 40 mins per neutrino energy on 512 
processors, and similar times for each of the two CC cross sections. 

The calculated cross sections for the NC- and CC-induced reactions are listed, respectively, in Tables III and IV 
for incident neutrino energies between 100 MeV and 1000 MeV. The columns labeled (1,C) and (1+2, C) [(1,PW) and 
(1+2, PW)] show results obtained by including fully interacting [plane-wave] NN final states, and one-body only or 
(one+two)-body terms in the weak current. These results are also displayed in Figs. 2-5. The two-body contributions 
are small, less than 10% over the whole energy range. Interaction effects in the final states are found to be even 
smaller, which suggests that realistic estimates for these cross sections on the deuteron (and possibly nuclei with 
A > 2) at high energies may be obtained by approximating the final nuclear states by plane waves, i.e. by employing 
the nucleon momentum distribution in the deuteron (or the spectral function in A > 2 nuclei). Lastly, in Figs. 2 and 4 
we also show the results obtained by including only the axial piece in the weak current. In this case, the interference 
response function R xy vanishes, and consequently the vi and V{ cross sections are the same. For CC-induced reactions, 
due to the charge-dependence of the NN final state (pp or nn) the neutrino-induced CC reactions has slightly larger 
cross section (a few %) even with only the axial piece in the weak current. We display the axial contribution in the 
anti-neutrino-induced CC reactions in Fig. 4. Axial contributions are larger than vector at low energy < (400-500) 
MeV, but become smaller than vector at higher energy. 

The sensitivity of the results on the model used for the two-body axial current (Set I or Set II) and NN potential 
(AV18 or CD-Bonn) is investigated, respectively, in Tables VI and VII. In both cases, the model dependence is found 
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TABLE III: Total cross sections in cm 2 for the NC-induced processes on the deuteron as function of the initial neutrino energy 
e, obtained with the AV18 potential and the inclusion of one-body (1) and (one+two)-body (1+2) terms in the weak current. 
Results corresponding to continuum final states (C) and plane wave final states (PW) are listed. 



e (MeV) 


H(1,C) 


n(l+ 


2,C) 


!/|(l,PW) 




1+2.PW) 


viihC) 


^(1+2,C) 






;i+2,pw) 


100 


2 


577 


-40) 


2 


657 


-40) 


2 


469( 


-40) 


2 


510 


-40) 


1 


604 


-40) 


1 


640 (-40) 


1 


607 


-40) 


1 


619( 


-40) 


150 


5 


720 


-40) 


5 


935 


-40) 


5 


626( 


-40) 


5 


752 


-40) 


3 


003 


-40) 


3 


075(-40) 


3 


096 


-40) 


3 


124( 


-40) 


200 


9 


435 


-40) 


9 


846 


-40) 


9 


384 ( 


-40) 


9 


650 


-40) 


4 


345 


-40) 


4 


460(-40) 


4 


526( 


-40) 


4 


576( 


-40) 


250 


1 


324 


-39) 


1 


389 


-39) 


1 


324( 


-39) 


1 


369 


-39) 


5 


531 


-40) 


5 


695(-40) 


5 


778 ( 


-40) 


5 


858( 


-40) 


300 


1 


683 


-39) 


1 


772 


-39) 


1 


689( 


-39) 


1 


754 


-39) 


6 


546 


-40) 


6 


762(-40) 


6 


842 ( 


-40) 


6 


962( 


-40) 


350 


2 


003 


-39) 


2 


116 


-39) 


2 


014( 


-39) 


2 


101 


-39) 


7 


420 


-40) 


7 


687(-40) 


7 


752 ( 


-40) 


7 


917( 


-40) 


400 


2 


279 


-39) 


2 


414 


-39) 


2 


295( 


-39) 


2 


403 


-39) 


8 


186 


-40) 


8 


504(-40) 


8 


545( 


-40) 


8 


760( 


-40) 


450 


2 


509 


-39) 


2 


664 


-39) 


2 


531( 


-39) 


2 


660 


-39) 


8 


856 


-40) 


9 


221 (-40) 


9 


255( 


-40) 


9 


520( 


-40) 


500 


2 


703 


-39) 


2 


874 


-39) 


2 


727( 


-39) 


2 


874 


-39) 


9 


503 


-40) 


9 


916(-40) 


9 


906 


-40) 


1 


023( 


-40) 


550 


2 


861 


-39) 


3 


046 


-39) 


2 


888( 


-39) 


3 


051 


-39) 


1 


010 


-39) 


1 


056(-39) 


1 


052( 


-39) 


1 


089( 


-39) 


600 


2 


989 


-39) 


3 


185 


-39) 


3 


019( 


-39) 


3 


196 


-39) 


1 


068 


-39) 


1 


118(-39) 


1 


110 


-39) 


1 


153( 


-39) 


650 


3 


093 


-39) 


3 


299 


-39) 


3 


125( 


-39) 


3 


315 


-39) 


1 


124 


-39) 


1 


178(-39) 


1 


166( 


-39) 


1 


214( 


-39) 


700 


3 


176 


-39) 


3 


390 


-39) 


3 


210( 


-39) 


3 


411 


-39) 


1 


178 


-39) 


1 


237(-39) 


1 


221( 


-39) 


1 


275( 


-39) 


750 


3 


243 


-39) 


3 


463 


-39) 


3 


278( 


-39) 


3 


489 


-39) 


1 


232 


-39) 


1 


295(-39) 


1 


275( 


-39) 


1 


333( 


-39) 


800 


3 


297 


-39) 


3 


522 


-39) 


3 


333( 


-39) 


3 


552 


-39) 


1 


284 


-39) 


1 


352(-39) 


1 


327( 


-39) 


1 


391( 


-39) 


850 


3 


340 


-39) 


3 


570 


-39) 


3 


377( 


-39) 


3 


603 


-39) 


1 


337 


-39) 


1 


408(-39) 


1 


379 ( 


-39) 


1 


448 ( 


-39) 


900 


3 


374 


-39) 


3 


608 


-39) 


3 


412( 


-39) 


3 


644 


-39) 


1 


388 


-39) 


1 


463(-39) 


1 


430( 


-39) 


1 


504 ( 


-39) 


950 


3 


403 


-39) 


3 


639 


-39) 


3 


440 ( 


-39) 


3 


678 


-39) 


1 


440 


-39) 


1 


518(-39) 


1 


481( 


-39) 


1 


559( 


-39) 


1000 


3 


425 


-39) 


3 


663 


-39) 


3 


461( 


-39) 


3 


704 


-39) 


1 


490 


-39) 


1 


572(-39) 


1 


530 


-39) 


1 


613( 


-39) 



TABLE IV: Same as in Table III, but for CC-induced processes. 



e (MeV) 


Fe(l,C) 


V 


e (l+2,C) 


V 


e(l,PW) 


z7 e (l+2,PW) 


I/e(l,C) 


V 


e(l + 2,C) 


V 


= (1,PW) 


f e (l+2,PW) 


100 


2 


567( 


-40) 


2 


606 


-40) 


2 


362 


-40) 


2.370( 


-40) 


6 


424 


-40) 


6 


631 


-40) 


5 


908 


-40) 


6 


023( 


-40) 


150 


4 


688( 


-40) 


4 


751 


-40) 


4 


487( 


-40) 


4.49K 


-40) 


1 


516( 


-39) 


1 


574( 


-39) 


1 


440 


-39) 


1 


477( 


-39) 


200 


6 


736( 


-40) 


6 


830( 


-40) 


6 


555( 


-40) 


6.568( 


-40) 


2 


605( 


-39) 


2 


719 


-39) 


2 


525 


-39) 


2 


603( 


-39) 


250 


8 


677( 


-40) 


8 


822( 


-40) 


8 


520( 


-40) 


8.567( 


-40) 


3 


775 ( 


-39) 


3 


958 


-39) 


3 


699 


-39) 


3 


833( 


-39) 


300 


1 


059( 


-39) 


1 


082( 


-39) 


1 


044( 


-39) 


1.056( 


-39) 


4 


928 


-39) 


5 


186 


-39) 


4 


854 


-39) 


5 


052( 


-39) 


350 


1 


254( 


-39) 


1 


286( 


-39) 


1 


239( 


-39) 


1.26K 


-39) 


5 


981( 


-39) 


6 


315 


-39) 


5 


923 


-39) 


6 


189( 


-39) 


400 


1 


455( 


-39) 


1 


499( 


-39) 


1 


441( 


-39) 


1.475( 


-39) 


6 


920( 


-39) 


7 


320 


-39) 


6 


876 


-39) 


7 


210( 


-39) 


450 


1 


663( 


-39) 


1 


722 ( 


-39) 


1 


650 


-39) 


1.698 


-39) 


7 


778 ( 


-39) 


8 


248 


-39) 


7 


704 


-39) 


8 


102( 


-39) 


500 


1 


879( 


-39) 


1 


952( 


-39) 


1 


865( 


-39) 


1.930( 


-39) 


8 


524( 


-39) 


9 


053 


-39) 


8 


410 


-39) 


8 


868( 


-39) 


550 


2 


100( 


-39) 


2 


189( 


-39) 


2 


087( 


-39) 


2.169 


-39) 


9 


064( 


-39) 


9 


636 


-39) 


9 


005 


-39) 


9 


519( 


-39) 


600 


2 


323( 


-39) 


2 


428( 


-39) 


2 


309( 


-39) 


2.410( 


-39) 


9 


556( 


-39) 


1 


017 


-38) 


9 


504 


-39) 


1 


007( 


-38) 


650 


2 


548( 


-39) 


2 


671 


-39) 


2 


537( 


-39) 


2.656( 


-39) 


9 


966 


-39) 


1 


062 


-38) 


9 


920 


-39) 


1 


053( 


-38) 


700 


2 


777( 


-39) 


2 


916( 


-39) 


2 


766 ( 


-39) 


2.905( 


-39) 


1 


031( 


-38) 


1 


098 


-38) 


1 


027 


-38) 


1 


091( 


-38) 


750 


3 


005( 


-39) 


3 


161 


-39) 


2 


995( 


-39) 


3.152( 


-39) 


1 


059( 


-38) 


1 


129 


-38) 


1 


055 


-38) 


1 


124( 


-38) 


800 


3 


232( 


-39) 


3 


403( 


-39) 


3 


223( 


-39) 


3.399( 


-39) 


1 


082( 


-38) 


1 


154 


-38) 


1 


079 


-38) 


1 


150( 


-38) 


850 


3 


456( 


-39) 


3 


645 ( 


-39) 


3 


448 ( 


-39) 


3.643( 


-39) 


1 


101 


-38) 


1 


176 


-38) 


1 


099 


-38) 


1 


173( 


-38) 


900 


3 


678( 


-39) 


3 


882( 


-39) 


3 


671( 


-39) 


3.885( 


-39) 


1 


118 


-38) 


1 


193 


-38) 


1 


116 


-38) 


1 


192( 


-38) 


950 


3 


896( 


-39) 


4 


116( 


-39) 


3 


890 


-39) 


4.122( 


-39) 


1 


131 


-38) 


1 


208 


-38) 


1 


129 


-38) 


1 


208( 


-38) 


1000 


4 


109( 


-39) 


4 


343 ( 


-39) 


4 


105 


-39) 


4.356( 


-39) 


1 


142( 


-38) 


1 


221 


-38) 


1 


141 


-38) 


1 


222( 


-38) 



to be negligible. The two-body vector currents are taken from the AV18, and therefore their short-range behavior is 
not consistent with the CD-Bonn interaction. This inconsistency, though, is of little numerical import. Furthermore, 
since interaction effects in the two-nucleon continuum appear to be negligible for neutrino energies > 100 MeV, 
the agreement between the calculated cross sections with the AV18 and CD-Bonn merely reflects the fact that the 
momentum distributions predicted by these two potential models are very close to each other for relative momenta 
< 400 MeV. 

In Figs. 6-14 we show the differential cross sections for the NC- and CC-induced reactions as function of the final 
lepton energy e' and scattering angle 8 at three incident neutrino energies, e = 100, 500, 900 MeV. The quasi-elastic 
peak is located at a final energy e' qc given by 

^ c = l + (2e/m) sin 2 0/2 ' (EU) 
where we have neglected the lepton mass in the case of the CC processes. Therefore as 9 changes from the backward 
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FIG. 2: (color online) Total cross sections for NC-induced processes on the deuteron, obtained with the AV18 potential and 
the inclusion of one-body (dashed line) and (one+two)-body (solid line) terms in the weak current. Also shown are the total 
cross sections obtained by retaining only the axial piece of the weak current. See text for explanation. 



TABLE V: (color on line) Total cross sections in cm 2 for the NC- and CC-induced processes on the deuteron obtained in 
Ref. [12] and in the present work at selected initial neutrino energies. Note that the values under the heading "this work" are 
slightly different than those reported in Table II for the reasons explained in the text. 



e (MeV) 


5 


50 


100 




Ref. [12] 


this work 


Ref. [12] 


this work 


Ref. [12] 


this work 


1/,-NC 
Pz-NC 
i/ e -CC 

T7 e -CC 


9.570(-44) 
9.364(-44) 
3.463(-43) 
2. 836 (-44) 


9.601 (-44) 
9. 403 (-44) 
3.440(-43) 
2. 842 (-44) 


5.944(-41) 
4.535(-41) 
1.376(-40) 
7.372(-41) 


5.942(-41) 
4.589(-41) 
1.367(-40) 
7.475(-41) 


2.711(-40) 
1.647(-40) 
6.836(-40) 
2.618(-40) 


2.703(-40) 
1.674(-40) 
6.735(-40) 
2.659(-40) 



to the forward hemisphere, the quasi-elastic peak moves to the right, i.e. towards higher and higher energies. Indeed, 
at forward angles it merges with the threshold peak due to the quasi-bound state. This latter peak is very 
pronounced at low e, but becomes more and more suppressed by the form factor ~ ( 1 So | jo(qr/2) \ d) as e increases. 

Finally, it is interesting to compare the results above with those obtained in a naive model, in which the deuteron 
is taken to consist of a free proton and neutron initially at rest. The lab-frame cross sections of the NC-induced 
processes on the nucleon, and of the CC-induced processes n(v e ,e~ )p and p(V e , e + )n in the limit in which the final 
electron/positron mass and proton-neutron mass difference are neglected, read [43]: 



da 
de'dfl 



NC/CC 



v j v 



G 2 m 2 
8tt 2 



- S(e'-e' c 



qeJ 



iNC/CC 



B 



NC/CC , ( s u ) 2 £<NC/CC 

m 4 



(5.2) 



where G—Gp or Gp cos8c for NC or CC, the — (+) sign in the second term is relative to the v (v) initiated reactions, 
e' qe has been defined in Eq. (5.1), and s — u = Ame — Q 2 with Q 2 = 4ee' sin 2 0/2. The structure functions A(Q 2 ), 
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FIG. 3: (color online) Total cross sections for NC-induced processes on the deuteron, obtained with the AV18 potential and 
the inclusion of (one+two)-body terms in the weak current. Also shown are the total cross sections obtained with plane-wave 
(PW) final states. 



TABLE VI: Total cross sections in cm 2 for the NC- and CC-induced processes on the deuteron at selected initial neutrino 
energies, obtained with the AV18 potential and the inclusion of one-body and set I or set II two-body terms in the weak 
current. 





vi- 


NC 


VI- 


NC 


v & - 


CC 


v E - 


CC 


e (MeV) 


set I 


set II 


set I 


set II 


set I 


set II 


set I 


set II 


5 

50 
100 


9.561 (-44) 
5.892(-41) 
2.657(-40) 


9.541 (-44) 
5.873(-41) 
2.652(-40) 


9.363(-44) 
4.546(-41) 
1.640(-40) 


9.344(-44) 
4.530(-41) 
1.636(-40) 


3.427(-43) 
1.348(-40) 
6.631(-40) 


3.421(-43) 
1.353 (-40) 
6. 621 (-40) 


2.831 (-44) 
7.403(-41) 
2.606(-40) 


2.826(-44) 
7.380(-41) 
2.600(-40) 



B(Q 2 ), and C(Q 2 ) for both NC and CC are given in terms of nucleon form factors in Appendix B. 

In the naive model, the v- and ^-deuteron NC cross sections are simply given by the sum of the corresponding 
proton and neutron (NC) cross sections, while the z/-deuteron (l^-deuteron) CC cross section is identified with the 
n(v e ,e~)p [p(V e , e + )n] cross section. The "model" differential cross sections as function of the final lepton scattering 
angle (after integrating out the energy-conserving (5-function) are illustrated in Fig. 15 at three incident energies 
(e = 100,500,900 MeV). The v and V cross sections are about the same at forward angles, for which Q 2 is small; at 
backward angles, as e and Q 2 increase, they both decrease due to the fall off in the form factors. However, this fall off 
is much more pronounced (orders of magnitude) for the V than for the v cross sections. (At low energy 100 MeV, the 
form factors do not change much with angle and the variation with angle in the differential cross section is mild, still 
it decreases more in the V than in v channel.) These features are reflected in the deuteron cross sections displayed in 
Figs. 6-14 (incidentally, in each panel of these figures the "model" cross sections would be represented by a ^-function 
located at e' qc , corresponding to the energy of the quasi-elastic peak). 

In order to illustrate nuclear correlation effects in the initial deuteron state, we compare the "model" v and V NC 
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FIG. 4: (color online) Same as in Fig. 2, but for CC-induced processes on the deuteron. 



TABLE VII: Total cross sections in cm 2 for the NC-induced processes on the deuteron at selected initial neutrino energies, 
obtained with the AV18 or CDB potentials and the inclusion of one-body terms (1) only and both one- and two-body terms 
(1+2) in the weak current. 









-NC 






v\ 


-NC 




e (MeV) 


AV18(1) 


CDB(l) 


AV18(l+2) 


CDB(l+2) 


AV18(1) 


CDB(l) 


AV18(l+2) 


CDB(l+2) 


50 
100 
500 
1000 


5.747(-41) 
2.577(-40) 
2. 703 (-39) 
3.425(-39) 


5.791(-40) 
2.597(-40) 
2.715(-39) 
3. 442 (-39) 


5.892(-41) 
2.657(-40) 
2. 874 (-39) 
3.663(-39) 


5.847(-40) 
2. 638 (-40) 
2. 858 (-39) 
3.659(-39) 


4.449(-41) 
1.604(-40) 
9.503(-40) 
1.490(-39) 


4.484(-40) 
1.617(-40) 
9.553(-40) 
1.496(-39) 


4.546(-41) 
1.640(-40) 
9.916(-40) 
1.572(-39) 


4.519(-40) 
1.633 (-40) 
9. 895 (-40) 
1.572(-39) 



cross sections with the plane-wave one-body results, shown in Fig. 16, for which we use the physical deuteron state, 
plane waves for the two-nucleon continuum states, and one-body currents. In both v and V NC reactions, inclusion of 
nuclear correlations in the initial state reduces the cross sections from the naive model. In fact, a similar reduction in 
cross section (due to nuclear correlations) at about nuclear density for uniform nuclear matter has been noticed, for 
example in Refs. [44, 45]. However, these correlations increase the ratio of v to V NC cross sections, as shown in the 
inset of Fig. 16. Similar effects are also found in the v and V CC reactions at low neutrino energy, as shown in Fig. 17. 
At higher energies, nuclear correlations hardly affect these cross sections, and the naive and realistic models are in 
better agreement with each other. The ratio of v to V CC cross sections is also increased due to nuclear correlations 
(see inset of Fig. 17). This fact may have interesting implications for long baseline neutrino experiments aimed at 
extracting CP violating signals from the detection of differences in the neutrino and antineutrino channels. 
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FIG. 5: (color online) Same as in Fig. 3, but for CC-induced processes on the deuteron. 

VI. CONCLUSIONS AND OUTLOOK 

In this work, we have studied inclusive neutrino scattering on the deuteron up to neutrino energies of 1 GeV, by 
using a realistic description of two-nucleon interactions and weak currents. Two-body terms in the latter increase 
the calculated cross sections with one-body currents by less than 10% over the whole energy region for both the NC- 
and CC-induced processes. Interaction effects in the two-nucleon continuum final state are found to be negligible for 
neutrino energies > 500 MeV. This suggests that fairly realistic estimates for these cross sections in light nuclei (and 
at relatively high neutrino energies) may be obtained in calculations based on the plane- wave impulse approximation. 
Even calculations in this limit, however, cannot be presently carried out, as they require knowledge of nuclear spectral 
functions over a wide range of missing momenta and energies — these are not yet available in light nuclei. Nuclear 
correlation effects in the initial deuteron state are found to be important. They reduce the v and, to a larger extent, 
V cross sections over the whole range of energies studied in this work, and therefore significantly increase the v to V 
cross-section ratio for both NC and CC reactions. In the present work the pion-production channels are not included. 
Experimentally they produce distinctive final states and give important contributions to total neutrino cross sections 
above pion-production threshold. It would be interesting to include these channels in future. 

It should be possible to use quantum Monte Carlo (QMC) methods [15] to study neutrino response functions, and 
associated sum rules, in light nuclei within the same (realistic) dynamical framework adopted here. Indeed, "exact" 
calculations of this type [46] led to a quantitatively accurate description of the quasi-elastic electromagnetic response 
functions measured in A = 3 and 4 nuclei. In particular, they showed that the charge-exchange character of the 
nucleon-nucleon interaction leads to shifts of longitudinal and transverse strength at higher excitation energies, thus 
providing a quenching of the response in the quasi-elastic peak region. This mechanism, however, is more than offset 
in the transverse channel by two-body currents, in particular those associated with pion exchange, and hence the 
response is enhanced over the entire quasi-elastic spectrum. It will be interesting to see the extent to which these 
considerations will remain valid in the weak sector probed in neutrino scattering, and possibly provide an explanation 
for the observed anomaly in the 12 C data. 
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FIG. 6: (color online) Differential cross section for NC-induced processes on the deuteron, obtained with the AV18 potential 
and the inclusion of one- and two-body terms in the nuclear weak current, as function of final lepton energy. The incident 
neutrino energy is 100 MeV. The final lepton angle is indicated in each panel. The (black) solid curve is for neutrino induced 
processes. The (red) dashed curve is for anti-neutrino induced processes. 
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Appendix A 



The cross section for CC processes at small incident neutrino energies in which the lepton mass cannot be neglected 
reads 



da 



G 2 k' 



F(Z,k') 



V j V 



voo -Roo + v zz R zz — v 0z R 0z + 

^xx-\-yy Rxx-\-yy ~F ^xy 



where the kinematical factors are given by 



k' 

v q = 2 e e' | 1 H — - cos ( 













Vzz = 


mf - 


h2ee' 


cos 6 


T 











+ ^K + 2c( e + e ')+<Z 2 ] , 



(Al) 

(A2) 
(A3) 
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FIG. 7: (color online) Same as Fig. 6, but the incident neutrino energy is 500 MeV. 
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mi is the final lepton mass, and the response functions are defined as in Eqs. (2.5)-(2.9). Note that 



(A4) 
(A5) 
(A6) 



e + e' = ^/2m 2 + (k + k') 2 + Q 2 , (A7) 
and the cross section above is easily shown to reduce to Eq. (2.3) in the limit m; = and Q 2 = 4ee'sin 2 9/2. 



Appendix B 

In this appendix, the structure functions entering the NC- and CC-induced processes on the nucleon are expressed 
in terms of (nucleon) form factors. In the NC case, they read: 



A NC = 4tj 
B NC = 



c nc = 



a,f n a {f n 1+ f n 2 ) , 

(^) 2 + (Ff) 2 + ,(Ff) 



N-N 



1 



(Bl) 
(B2) 
(B3) 
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FIG. 8: (color online) Same as Fig. 6, but the incident neutrino energy is 900 MeV. 



and r\ = Q 2 / (Am 2 ). The nucleon form factors if and F A for N = p or n are defined as 

2 F v J n = (1 - 4 sin 2 9 W ) Ff /n - F? /p , (B4) 
2F P A /n = T Ga, (B5) 

where the proton and neutron electromagnetic form factors are, respectively, Ff = (Ff+F^)/2 and F™ = (Ff — F^)/2 

with F^ V defined in Eqs. (3.43)-(3.44), and the axial form factor Ga (with - for p and + for n) as defined in Eq. (3.12). 
In the limit in which the final lepton mass and proton-neutron mass difference are both neglected, the relations for 
the A, B and C structure functions remain valid for the CC case, provided 

if — ► FY , F N A -^G A . (B6) 
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FIG. 10: Same as Fig. 9, but the incident neutrino energy is 500 MeV. 
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FIG. 11: Same as Fig. 9, but the incident neutrino energy is 900 MeV. 




FIG. 12: Differential cross section for electron neutrino induced CC processes on the deuteron, obtained with the AV18 potential 
and the inclusion of one- and two-body terms in the nuclear weak current, as function of final lepton energy. The incident 
neutrino energy is 100 MeV. The final lepton angle is indicated in each panel. 




FIG. 13: Same as Fig. 12, but the incident neutrino energy is 500 MeV. 
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FIG. 14: Same as Fig. 12, but the incident neutrino energy is 900 MeV. 
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FIG. 15: (color on line) The "model" NC (left panel) and CC (right panel) differential cross sections for neutrino (solid lines) 
and antineutrino (dashed lines) energies of 100, 500, and 900 MeV, as functions of the final lepton scattering angle. 
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FIG. 16: (color online) The "model" (P+N) NC cross sections for neutrino and antineutrino are compared with plane-wave 
one-body (PW 1-body) results, see text for explanation. Inset: ratio of neutrino NC versus antineutrino NC cross section. 




FIG. 17: (color online) Same as Fig. 16, but for CC cross sections. 



